Reconstruction of Hamiltonians from given time evolutions 



o 

(N 

IX, 

l> 
(N 



> 
oo 
in 

O 
O 



X 



Reconstruction of Hamiltonians from given time evolutions 

J. Bernatska 1 '^ and A. Messina 2 'E3 

^National University of 'Kiev-Mohyla Academy', 2 Skovorody Str., 04070 Kiev, 
UkrainW\ 

2 ' Dipartimento di Fisica, Universitd di Palermo, Via Archirafi 36, 90123 Palermo, 
Italy 

(Dated: 1 March 2011) 

In this paper we propose a systematic method to solve the inverse dynamical problem for a finite-level quantum 
system governed by the von Neumann equation: to find a class of Hamiltonians reproducing a prescribed time 
evolution of a pure or mixed state of the system. Our approach exploits the equivalence between the action 
of the group of evolution operators over the state space and the adjoint action of the unitary group over 
Hermitian matrices. The method is illustrated by two examples involving a pure and a mixed state. 
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I. INTRODUCTION 

The problem of reconstructing a Hamiltonian or a La- 
grangian from observations on a physical system is a fun- 
damental issue in physics. Its history can be traced back 
to the time of Kepler and Newton when the law of grav- 
itation was discovered from observations of the plane- 
tary motion. The reconstruction of an effective potential 
from the data of particle scattering also belongs to this 
kind of problems. In the context of quantum control the- 
ory reconstructing Hamiltonians on demand represents 
a strategic intermediate step toward the realization of 
noise-protected control protocols governing the dynam- 
ics of a given, often complex and multipartite, quantum 
system. 

To drive the time evolution of a quantum system at 
will is at the heart of the quantum control theory. Fruit- 
ful applications of this theory in physical chemistry al- 
low to achieve the optical control over ultrafast molecular 
processes^, and to design feedback control for quantum 
state preparations in quantum optics^. A specific task of 
the quantum control theory is to contrive schemes that 
allow to steer a quantum system to a desired state start- 
ing from a given initial one. More ambitiously one could 
even require a possibility to control an evolution of the 
system under scrutiny at any instant of time. In other 
words, one could require to design the control in such a 
way that the system evolves along a prefixed path in its 
state space. Being able to solve this problem, a dream 
since the early days of quantum mechanics, would help 
to highlight dynamical behaviors of microscopic systems, 
and would provide additional tools and ideas applica- 
ble to quantum computation. In a recently published 
paper— the tomographic method has been used to show 
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that a specially modeled time-dependent state of two 
qutrits, violating Bell's inequalities at t = 0, loses this 
non-classical feature at successive times. In this paper 
we introduce a systematic method to solve the follow- 
ing inverse problem of the unitary quantum dynamics: 
given the time evolution of a system is explicitly speci- 
fied by means of a density matrix p, to construct a class 
of generally time-dependent Hamiltonian models whose 
associated von Neumann equations are satisfied by this p, 
describing a pure or a mixed quantum state. 

The paper is organized as follows. In Section 2 math- 
ematical notations are introduced. In Section 3 we ex- 
plain our mathematical procedure: how to reconstruct 
the Hamiltonian reproducing the prescribed time evolu- 
tion, if this evolution obeys the von Neumann equation. 
Every time evolution is realized by a class of Hamiltonian 
models. We illustrate our method in Section 4 where it 
is applied to time evolutions of both a pure state and a 
mixed one. 



II. MATHEMATICAL NOTATION 

To establish our notation we start with brief recalling 
the mathematical grounds of quantum mechanics after 
J. von Neumann^— . A quantum system is generally de- 
scribed by three basic ingredients: states, observables, 
and dynamics (the law of time evolution) . Very often the 
state space is identified with a complex Hilbert space ri. 
In what follows we consider finite-level systems, so the 
Hilbert space is finite dimensional. After Dirac we de- 
note its vectors by |$), \*&), . . . , and an inner product 
by ($|*). 

But the Hilbert space represents only pure states. 
Every vector |$) S H gives rise to the ray (the one- 
dimensional subspace) {A|<!>) | A <E C}, associated with 
one pure state, say ip. In other words, a pure state is 
identified with an equivalence class of vectors of length 1 
in "H, and two vectors represent the same pure state if 
they differ only by a phase factor. 
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If we want to include mixed states in our pattern, it 
is better to represent the pure state ip as a projection 
operator P$ onto the ray {A|$) | Ae C}. In Dirac's no- 
tation the operator Pj> is expressed as |$)(<I>|. Every 
such projector Pj>, and so every pure state <p, is a self- 
adjoint linear transformation on the Hilbert space such 
that P§ = P$ . The collection of all pure states form a 
projective space. A 'statistical mix' of pure states 

P = J2PkP^ k , 

k 

where ^ pu < 1 and ^2 k pi~ = l, represents a mixed state. 
The operator p is not a projector, but a self-adjoint pos- 
itive definite linear transformation. Therefore the state 
space of a quantum system can be identified with the 
subspace of positive definite matrices in the space 21(H) 
of self-adjoint linear transformations on H. In order to 
identify a quantum state with its density matrix we re- 
strict the state space to matrices of trace 1. 

We introduce a scalar product in 21(H), and the dual 
space to 21(H) as follows: 

(p, A) = Tr pA, p€ 21(H), A e 2T(H). 

The dual space 21* (H) represents the associative algebra 
of observables^ of a quantum system. The expectation 
value of an observable A with respect to a state p is 
defined by the formula 

(A) p = (p,A)=TrpA. 

Then for a pure state ip = Pj> we have 

{A) v = TrP*A={$\A$). 

The last ingredient of a quantum system is its dy- 
namics, which is governed by the well-known Schrodinger 
equation 



ih 



dt 



H\9>) 



m en, 



if we deal with the pure state associated with the vec- 
tor |$); or in general by the von Neumann equation 

ih^ = [H,p], pe 21(H), 

where [■, •] denotes a commutator. The self-adjoint oper- 
ator Hg21*(H) is called a Hamiltonian. 

According to Stone's theorem the time evolution gen- 
erated by a timc-indcpcndcnt self-adjoint operator H is 
realized by unitary operators {£/*}, forming a continu- 
ous one-parameter group in t, that is UtU s = Ut+ s for 
all t, s £ M. and lim^o Ut—>I- Then the Schrodinger 
equation with an initial condition at to is solved by the 
formula |$(t)) = U t -t |3>(*o)) ■ In the case of a time- 
dependent H time evolution is also realized by a con- 
tinuous group of unitary operators^ {Ut,t } such that 



Ut a ,tiUt lt t = U t3t t , Ut ,t = I; and the Schrodinger equa- 
tion is solved by =E/ Mo |$(i )>. 

Hereinafter wc deal with a more general pattern, where 
the state space is identified with the subspace V of pos- 
itive definite matrices of trace 1 in the space 21(H), and 
dynamics is governed by the von Neumann equation. 
Time evolution is realized by the group of unitary op- 
erators {U tt t } as follows 



p(t) = u t , toP (t )u^ . 



(1) 



The Hamiltonian generating this evolution relates to the 
evolution operator £/t,t according to the formula 



H(t) = ih^ 



u; 



(2) 



which can be obtained by substituting the evolution of 
| $) , or p, into the corresponding differential equation. 

The mentioned well-known facts of quantum mechan- 
ics are naturally interpreted from the position of the Lie 
group and algebra theory. Evidently, time evolution of a 
quantum system is realized by the adjoint action (fTJ of 
a unitary group. This gives us a hint how to treat the 
problem. Suppose we deal with a n-level system, then 
the state space can be identified with the positive definite 
subspace V of the affine linear space isu(n) © (l/n)I„. In 
our context the latter space can be replaced by the alge- 
bra isu(n), what we show in the next section. Therefore, 
from the mathematical point of view we deal with the 
adjoint action of a Lie group over its Lie algebra. 

Moreover, it is easy to check that Tr p k remains invari- 
ant under time evolution for any k 1: 



Tr p k (t) = Tr / (0) = const, 



(3) 



where only first n — 1 equations are functionally inde- 
pendent for a n-level system. In other words, the adjoint 
action ([1]) divides the state space into non-intersecting 
orbits, and every orbit is determined by the set of equa- 
tions ([3]) . This clarifies the geometry of time evolution of 
a closed quantum system. Every quantum state evolves 
within a single orbit, fixed by an initial state. It means: 
only states within this orbit are reachable from the initial 
state. 



III. SCHEME OF RECONSTRUCTION FOR n-LEVEL 
SYSTEMS 

We aim at solving the problem of reconstructing a 
Hamiltonian model from a given time evolution satisfy- 
ing the von Neumann equation. Actually we reconstruct 
an evolution operator, which is considered as an element^ 
of a unitary group. For this purpose we use a generalized 
stereographic parametrizatio n 10 ^ 1 , applicable to states 
with finite number of components, say n. Therefore in 
what follows we deal with matrices of size n x n, and 
mark them by the hat sign ~. 
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A. An orbit structure of the state space 

First of all, we show an equivalence between the de- 
scribed above physical and mathematical patterns. Let 
us consider the algebra isu(n) consisting of Hcrmitian 
traceless matrices. We denote these matrices by /t, keep- 
ing p for matrices of trace 1. The adjoint action of the 
unitary group SU(n) over isu(n) is naturally defined by 

ji = u(ii n u~ , p, £ isu(n) , MGSU(n). (4) 

All [i obtained from the initial point Am when u runs the 
group SU(n) form an orbit. This orbit can be determined 
by the set of equations 

Tr fi k = Tr A m = const, k = l, ... n-l. (5) 

The whole algebra isu(n) is devided into non-intersecting 
orbits by the adjoint action of SU(n). 

It is easy to check that the adjoint action over isu(n) 
in (01 and the time evolution over the state space V in 
(fTj) are realized by the same unitary group. It means that 
for a n-level quantum system the group of unitary oper- 
ators {Utt } i s SU(n). Also the algebra isu(n) has the 
same orbit structure as the affine space isu(n) © (l/n)I n 
containig the state space of the system. Indeed, every or- 
bit in isu(n) determined by ([5]) corresponds to a unique 
orbit ^ in isu(n) ffi (l/n)I n . Therefore instead of the 
affine space we consider the Lie algebra isu(n). Also we 
shift the state space V into the corresponding subspace 
M. Cisu(ro) by subtracting the matrix (l/n)I„. 

Now we can use the knowledge of the orbit structure of 
a Lie algebra under the adjoint action of its Lie group in 
order to describe the structure of the state space. Above 
we called this structure the geometry of time evolution. 
Every orbit in isu(n) is specified by its initial point taken 
from a fundamental Weyl chamber or its walls. The Weyl 
chamber is a subset of the Cartan subalgebra f) of isu(n), 
which is defined as the maximal Abelian subalgebra of 
a semisimplc Lie algebra. For isu(n) it consists of real 
diagonal matrices of zero trace. Thus a canonical form of 
A € isu(n) can serve as an initial point for the orbit where 
A is located. But there is an ambiguity upon the order 
of diagonal entries. This ambiguity is removed by intro- 
ducing the Weyl group, which performs a permutation of 
the diagonal entries, and whose fundamental domain is 
called a fundamental Weyl chamber. 

Here we choose the positive Weyl chamber as a funda- 
mental one. Its interior is the open domain 

C = {A e t) | (A, > 0, Va G A+}, 

where A + denotes the set of positive simple roots. A wall 
of the positive Weyl chamber is the set 

r„ = {Aef) I (P,a) =o}, aeA+ 

or any intersection of such sets. 

We omit the accurate definition of simple roots, instead 
we represent them in the standard for the algebra isu(n) 



form. The simple roots belong to the dual subalgebra rj* 
and form a basis in it. Because f)* is equivalent to f) we 
can use the same representation. The set A + consists of 
n—1 diagonal matrices {&k = Ef e i c —Ek+i,k+i}i where the 
matrix Ejk has the only nonzero entry in row j column k 
and this entry equals 1. It is easy to construct the dual 
basis {7 (fe) } C f): 

k 

Then every point of the positive Weyl chamber of zsu(n) 
can be represented in the form 

n-1 

Ain = $> fe 7 (fe) , b k em+. 
fc=i 

If we allow the coefficients {bk} to equal zero, we obtain 
the closed positive Weyl chamber containig all its walls. 
Every wall spans a nontrivial subset of {7^}, and is 
represented by a hyperplane bounding the chamber. 

Points from the interior of a fundamental Weyl cham- 
ber in isu(n) give rise to generic orbits of dimension 
n 2 — n. Every point of a wall of the Weyl chamber gives 
rise to a degenerate orbit of lower dimension. In the case 
of a onc-dimcnsional wall one obtains a maximal degen- 
erate orbit of dimension 2(n— 1). The latter represents 
a set of pure states if being shifted back into the space 
isu(n) © (l/n)I n remains a projective space. 

In particular, the points bi^ 1 ^ and 6„_i7 ( -™~ 1 -' give 
rise to maximal degenerate orbits in Ai if b\ ^ 1, 
b n -i ^5 l/( n ~ !)■ But only one of them, namely 7*- 1 -*, 
generates the set of pure states. Positive definiteness of 
the state space V implies that acceptable values of {6^} 
satisfy the inequalities 

n— 1 n—1 

b k - ^ kb k + 1 > 0, j = l,...,n. 

k—j k—1 

The last inequality (j = n) appears to be sufficient if one 
takes into account non-negativity of {bk}. 

Obviously, the elements of the basis {7^}, except 7*- 1 ', 
do not belong to the subspace Ai. We are able to con- 
struct a basis in f) fulfiling this condition. This new basis 
consists of traceless real diagonal matrices {-7r( fc )} such 
that (tt^) 2 — 7r( fe ) = 0. It is easy to write these matrices 
cxplicitcly: 

7r (fc) = E kk - (l/n)in. 
The basis {nk} consists of projectors, and the expansion 

n-1 

Am = Ck%ih) ( 6 ) 
fe=l 

immediately displays a dimension of the orbit passing 
through Ain- The dimension is (n — l)(m+ 1), where m 
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denotes a number of nonzero summands. Simply speak- 
ing the degeneracy of an orbit becomes evident from co- 
incidence of diagonal entries of pi n . For example, if all 
but one entries are equal, then the orbit passing through 
/tin is maximal denerate. 

Therefore, every evolution of a closed n-level quantum 
system is located within a certain orbit of the adjoint 
action of the group SU(n). The orbit is fixed by an initial 
state, say p(t ). It means that only the states within the 
fixed orbit are reachable from p(to) because orbits do not 
intersect. All orbits together form the state space V. 

Turning back to our problem, suppose we are given 
a sample of evolution p as a function in £, and we aim 
at constructing an operator Ut.t £ SU(n) realizing the 
given evolution. First we should identify the orbit where 
p lives. We propose to identify it by the canonical form 
of p, say pj n . The matrix p; n is constant and character- 
izes the whole orbit. In particular, it serves as a canon- 
ical form for p{to). The constancy of p m follows from 
the set of equations ([3]) determining an orbit. Indeed, 
fixed values of Tr p k guarentee that all poins of an or- 
bit have the same characteristic polynomial, and so have 
the same eigenvalues. In order to detect a degeneracy of 
the orbit, one should apply the expansion (J6j> to the ma- 
trix /5; n — (l/n)I„, that displays a dimension of the space 
where the given state p evolves. This determines a struc- 
ture of the evolution operator Ut,t an d some auxiliary 
matrices in the proposed method. 

The next step of our scheme consists in constructing 
a generic form of the unitary matrix u <= SU(n) realizing 
motion within a chosen orbit. We call this matrix an 'evo- 
lution' one, indicating by quotes the difference between 
a motion and a time evolution realised by an evolution 
operator Ut,t a - 



B. Stereographic parametrization and 'evolution' 
operators 

As shown above, motion of a n-level quantum system 
is realized within an orbit by elements of the unitary 
group SU(n). We call these elements 'evolution' matri- 
ces when they concerns motion within orbits, and evolu- 
tion operators when they generate time evolution. There 
is a connection between an orbit and its 'evolution' ma- 
trix. An 'evolution' matrix can be written in terms of 
the complex coordinates parameterizing an orbit. In the 
paper we implement this by the method of generalized 
stereographic parametrization. 

First wc briefly recall the idea of the generalized stere- 
ographic parametrizatio n 10 ' 1 1 for a compact Lie group G. 
It consists in comparing the Iwasawa and Gauss decom- 
positions of the complexified group G c whose maximal 
compact subgroup is G. Let T be the maximal torus of 
G, and T c be the maximal torus of G c . We also define 
the following subgroups of G c : the maximal real Abelian 
subgroup A, a nilpotent subgroup R, and another nilpo- 
tent subgroup Z = R* . Note that T c = T A. The Iwasawa 



and Gauss decompositions of G , respectively, are 

G c = GAR, G c = ZT C R. 

If we define a generic adjoint orbit by means of coset 
space we obtain the equality line 

G _G C ZT C R 
= T ~ ~P ~ T C R ~ ' 

where P = TAR is the minimal parabolic subgroup of G c . 
Therefore, we can parameterize the orbit O in terms of 
canonical coordinates of Z, we denote them {zfc}. 

In general an orbit of the adjoint action of group G 
is defined as the coset space G/G Pia , where G Pin denotes 
the stability subgroup of the point pi n identifying the 
orbit. A generic orbit arises when a stability subgroup 
coincides with the maximal torus T. A larger stability 
subgroup gives rise to a degenerate orbit, then a nonmin- 
imal parabolic subgroup P = G Pin AR should be taken. 
In this case the nilpotent subgroup Z is generated by a 
smaller number of root vectors, and is parameterized by 
a smaller number of canonical coordinates. 

Here wc deal with the group SU(n). Its complexifi- 
cation having SU(n) as a maximal compact subgroup is 
SL(n,C). All mentioned subgroups can be represented 
in matrix form. The Abelian subgroup A consists of di- 
agonal matrices a = diag(c?i, d,2, . . . , d n ) with J| fe dk = 1. 
Let R be a subgroup of upper triangular matrices r , then 
Z is a subgroup of lower triangular matrices z. 

For a matrix zG Z, assigned to a point of an orbit, we 
obtain an 'evolution' matrix u £ SU(n) from the Iwasawa 
decomposition of z: 

z~uaf. (7) 

At the same time we aim at expressing the matrix u 
in terms of the canonical coordinates {zfc} apperaing as 
entries of the matrix z. The easiest way to achieve this 
goal is to solve the equation 

z z = r a u uar = r a r. (8J 

for all entries of a and f. Having the matrices a and 
f expressed by means of {zk}, one easily calculates the 
matrix u from the equality ([7]). 

As a result, we obtain a generic form for the ma- 
trix u, parameterized by the complex canonical coordi- 
nates {zk}. This matrix realizes the 'evolution' starting 
from zero-point {zk = 0} within every orbit of the same 
type. 

Remark 1. It should be noticed, that the obtained 'evo- 
lution' matrix u is not a unique solution of the problem. 
One can multiply it by any matrix v € G Pin , that gives 
another 'evolution' within the same orbit. 



C. Pure states 

Here we consider how to apply the above approach to 
the particular case of pure states. Usually such states are 
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described by vectors of the Hilbert space H. Let {|A;)} 
denote a set of basis quantum states, then a pure state 
can be expressed in the following form: 



W)) 



k 



Ck(t)\k), 



where the coefficients {ck} satisfy the normalization con- 
dition \ c k(t)\ 2 = 1- One can reduce this expression to 
the form of a projective vector by introducing the com- 
plex coordinates Zk=Ck/c\ for k = 2, n which pa- 
rameterize the maximal degenerate orbit where the state 
lives, and z £ Z has the following form 



/ 1 

Z 2 1 

z 3 1 



\z„ 

Let f and a have the form 












(\ 



ri2 r 13 
1 r 23 
1 



\ 

(a 2 

a 3 /a 2 
04/03 



V 











r 2n 




r 3n 




1 J 




\ 








l/a, 


J 



After solving JSJ we obtain the expressions 



= 1 



El* 



Zk/al, k ^ 2, 



-ZjZ k /a 2 J+1 , + 



which provide a parameterization of the matrices a and f. 

The matrix u calculated by (J7J rcallizcs evolution of the 
pure state \ifj) starting from the initial state \tp(to)) = |1). 
In the case of other initial state \ipo)i one should multi- 
ply u by the matrix Uq 1 such that uq realizes the evo- 
lution form |1) to IV'o)- Having the evolution operator 
Ut.to == "'"(7 1 : one computes the corresponding Hamilto- 
nian by the formula 



6 = ihS*L 
dt 



U ttt0 = ih-u 



which does not depend on the initial state. 

The obtained evolution operator Ut,t realizes exactly 
the prescribed evolution of but is not a unique so- 
lution of the problem in question. According to Re- 
mark [T] any evolution operator Ut,t = uvv,q , where 
v £ U(l) x SU(n — 1), realizes the same evolution of 
If v is a function in t, the corresponding Hamiltonian 
changes into 



m = h 



„ dv „ _ 
ihu—v 
dt 



(9) 



which is a gauge transformation of H. 



D. Mixed states 

Considering mixed states we represent them by density 
matrices, which are Hcrmitian positive definite matrices 
of unit trace. Because an orbit is specified by a diagonal 
matrix, one should diagonalize the given matrix p: 

p = D Pin D-\ 

where pi n is the canonical form of p. As proven at the 
end of section ITlIAi an orbit structure of the state space 
implies that the matrix pi n is constant and serves as the 
canonical form of p(to). Preserving the Hcrmitian char- 
acter of p, the matrix D should be a unitary operator, 
so it can serve as the desired evolution operator realizing 
the given evolution p. Here we show how to obtain such 
operator by means of the proposed method. 

We start with constructing a generic element the group 
SU(n) according to the scheme of subsection IIII Bl from 
the matrix 



/ 1 

Z12 

Zl3 




f 

z 23 



\Zln 



z 2n z 3n 








1/ 



In order to find relations between the entries of p and the 
canonical coordinates {zjk} we solve the equation 

p = upinu' 1 

for {zjk}- The relations allow to introduce the prescribed 
time evolution into the matrix u. The obtained time 
dependent matrix u realizes the evolution starting from 
pin, If p(to) does not coincide with its canonical form p 
the matrix u should be multiplied by the matrix u\ 



1 



such that /3(io) = &o/3i> 



Finally, the given evolution 



from p(to) to p(t) is realized by the operator Ut,t Q = i 
Evidently, the Hamiltonian does not depend on uq: 



tt .f- dU t,to T y-i 1 
H=ih — —U+t = ih — u . 
dt t,io dt 



We obtain the Hamiltonian model guaranteeing that p 
evolves according to the prescribed law. Again this model 
is not unique. Any evolution operator U t 



t,t 



'Uq with 

the matrix v belonging to the stability subgroup of the 
orbit where the state p lives also realizes the prescribed 
evolution of p. In this case v performs a gauge trans- 
formation of the model, where the Hamiltonian changes 
according to the law 



IV. ILLUSTRATION OF THE METHOD 

To illustrate the method proposed in the previous sec- 
tion we consider two examples: a pure state and a mixed 
one. In the both cases we calculate evolution operators 
and the corresponding Hamiltonians, and show how to 
construct other possible evolution operators. 
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A. The case of a pure state 



gives 



As a time dependent normalized pure state we choose 
the following one^ 



\m) = 



i 



■{|l,-l)+cosw<|0,0) + 
V2 + cos^ Lot 

,1)}, (10) 



i7r/6| 



-1, 



describing a time evolution of a pair of interacting spins 1: 
S\ and S 2 - The state \mi,m2) is a simultaneous eigen- 
state of Siz and S 2z , that is Si Z \mi,m 2 ) =mi\mi,m 2 ), 
where m; = — 1, 0, 1 for 7 = 1, 2. 

It is worth noting that the evolution (|10p of \ip(0)) 
exhibits the conservation of m = (mi + m 2 ) = 0, which 
is an eigenvalue of the observable S z = Si z + S 2z . As 
a result, \ip), which is indeed a vector from the nine- 
dimensional Hilbert space T-Li 2 of the two qutrits, lives in 
a three-dimensional subspace 77o CH.ii of S x pertaining 
to to = 0. Thus, we need not to take into consideration 
the whole Hilbert space 7^12 ■ It is enough to describe the 
evolution of \tp) by a reduced Hamiltonian model living 
in Hq. 

Another interesting observation is the following. Fo- 
cusing on the mathematical nature of the problem 
we could replace the specific evolution (JTUJ) of two 
qutrits by the evolution of a single qutrit J such that 
J Z \M)=M\M) with M = 0, ±1: 



\m) 



1 



v2 + cos 2 Lot 



+coswt|0) +e 47r / 6 |-l)} 



(11) 

Both evolutions: ((TU)) and (fTTj) can be described by the 
same Hamiltonian model constructed below. 

For the sake of generality we consider a state of the 
form 

\m) = ci(t)\l,-l) + c 2 (t)\0,0) + c 3 (t)\-l,l) (12) 



with \ Cl (t)\ 2 + \c 2 (t)\ 2 + \c 3 {t)\ 2 = 1. While CCD is a pure 
state, its evolution is located within a maximum degen- 
erate orbit SU(3)/(U(1) x SU(2)), which is the projective 
space CP 2 . Reduction to the form of a projective vector 



m)) = |l,-l) + «a(i)|0,0) + z 3 (*)|-l,l>, 

where 2:2 = C2 j c\ , 23 = C3 j c\ . On the other hand, this 
projective vector can be rearranged into the matrix form: 




From we obtain the following: 



*2 



1 



'12 



T23 



ri3 



Then an evolution matrix u looks like 



a 2 0,2113 
22 03 



ci 2 

Z3 

v C12 0,2(13 



a 2 _ 

- -2 



(13) 



This matrix acting on the vector e = (1,0,0) gives 
z= (l/a2, ^2/02, 23/(12); where 02 serves as its norm. 
For an arbitrary initial vector zq we need to find the 
matrix ii such that UQe~z . Then 

Z = UUq 1 Zo, 

and U t = uUq 1 is the required evolution operator. 

Now we come back to the state (jlOjl . In this case, 
z 2 (t) = cos Lot, Za(t) = e I7r / 6 , and a 2 = 2 + cos 2 ut, a\ = 2. 
Then 

/ 



V2+cos 2 ujt v /2(2+cos 2 UJt) 
COS U)t y/2 



\/2+cos- ojI 

in/6 



\/2+cos 2 



y V2+cos 2 ojt ^/2(2+cos 2 wt) 

/ 




1 

^2 



u 



1 


-1 






1 


v'2 


vs 


y/S 




_ e ™/6 



V V3 

and the evolution matrix is 



V6 



V2 



1 

v^2 



I 2 + cos ut+-^3(2 + cos 2 Ljt ) 

2-y/3(2+cos 2 ut) 

cos ut — 1 
2-y/3(2+cos 2 wt) 
2+cos ut--^/3(2+cos 2 ut) 
\ 2-^/3(2 + cos 2 ut) 



t/6 



^/3(2+cos 2 ut) 

2 + cos uit 
-y/3(2+cos 2 ut) 
1 — cos ut rt^/S 



^/3(2+cos 2 ut) 



2+cos ut — ^/ 3 (2+cos 2 ut) — z7r/6\ 
.-. /TTTTT", t 77 ^ 



2^/3(2+cos 2 ut) 

cos ut — 1 
2^/3(2+cos 2 ut) 



-Z7T/6 



2+cos ut+-\/3(2+cos 2 ut) 
2^/3(2 + cos 2 ut) 



The Hamiltonian calculated by @ is the following and coincides with the result reported in^. 

ihio sin tot ( ^ ^ ^ \ As mentioned in Remark[I] the obtained evolution ma- 

7? = - — 1—1 — e l7r / 6 i (14) trix J7 t is not a unique solution of the problem. As well 

+ cos Lot \ q e -i7r/6 / the problem is solved by any matrix U t = uvuQ 1 with 7} 
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of the following form 



,iy(t)/H 








a(t) (3{t) 
-/?•(«) a*(t) / 



where 7 is an arbitrary real-valued function in t, and 
a, j3 are complex- valued functions in t such that 
|a(i)| 2 + |/3(i)| 2 = l for all t. 

It is easy to check that the reduced to Ho^ operators 
Su+S 2z , Si z S 2z , StS 2 , and S^Sf, where Sf=S ix ± 
iSiy, i = 1, 2, can be cast in the following matrix form: 



st + si 



St S 2 — 




Si si = 



Si S 2 — 




(15) 



Exploiting the expressions (|15j) a Hamiltonian model 
H gcn acting in T-Li 2 , whose reduction to Ho is given by 
(fill) , is written in the following form 



H gcn — htt Q (Sf + Sf ) + 
ih\(t)({e- l iStS2 + SiS+)SfSi 

SfSKStSz +e^S~St)' 



(16) 



(17) 



where the coupling constant A has the form 

1 u> sinwi 
' ' ~ ~2 ' 2 + cos 2 ujt 

and MIq is an arbitrary value of energy. Equation f| 16[) 
may be interpreted as describing two interacting qutrits 
Si and S 2 subjected to a magnetic field. The Hamil- 
tonian H gcn commutes with S z at any t and guarantees 
that the value of m— (mi +m 2 ) remains constant, thus 
the evolution of \ip(0)} belongs to a maximal degenerate 
orbit of SU(3) over the state space. 



B. The case of mixed state 

To illustrate our method in the case of a mixed state, 
let us consider a two-level atom coupled to a bosonic 
mode. Atomic dynamical variables can be represented 
in terms of the Pauli matrices a x , <j v , a z acting upon a 
two-dimensional Hilbcrt space where a basis is {|+), |— )} 
with er z |±) = ±|±)- The bosonic mode is represented, as 
usual, by an annihilation a and a creation operators 
such that [a, a'] = 1, acting upon an infinite-dimensional 
Hilbert space with the basis {|n),n = 0,1,2,...} such 
that a^a\n) = n\n). Suppose that at t = the composed 
system is in the mixed state 

p(0) = (3/4)|0, +}(0, +| + (1/4)|1, (18) 



where \n, ±) is a simultaneous eigenstatc of a^a and a z 
pertaining to their eigenvalues n and ±1 respectively. As- 
sume that the atom-mode interaction generates the fol- 
lowing evolution: 

p(t) = (3/4-(l/2)cos 2 5 <)|0,+)(0,+| + 

(l/4+ (1/2) co S V)|l,-> <1,-| + 

(1/4) sin2 5 t(|0, +)(1, -I + |1, -)(0, +|). (19) 

It is suitable to represent p in a matrix form as follow: 

^3/4- (1/2) cos 2 5 t (l/4)sin2.g£ \ 
P[ >- \ (l/4)sin2<7i 1/4+ (1/2) cos 2 3^ ' 1 ' 

which may be generalized by the following notation for 
entries of p: 



P 



P11 P12 



K Pl2 P22 

Diagonalization of p leads to 

Ai 
A 2 



Ai,2 = |( P11 + p 22 ± \/(pu ~ P22) 2 +4|pi 2 | 



The group acting over the state (f20|) is SU(2), and its 
generic element obtained by the explained above method 
of stereographic parametrization has the form 



1 -z 



y/l + \z\ 2 V 1 

Solving the equation u/5i n u _1 = p for z, we come to 



2p 



12 



P11 - P22 ± y/ipn ~ P22) 2 + 4|/9i 2 | 2 
In the case of state (|20p we obtain 



'1/4 
pin "I 3/4 



and z — — tangt. Then a possible evolution operator and 
its correspondent Hamiltonian acquire the forms 

t fuxgt sin ff A ^ = / l\ 
\~svngt cos gt J ' \ _1 °/ 

The given evolution ([20]) is obtained by 

p(t) = u t mu t - 1 - 

According to Remark [TJ there exist other solutions of 
the problem. Evolution operators Ut = uv such that 

e ia(t)/h Q 
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where a and (3 are arbitrary real- valued functions in t, 
also realizes the same evolution (f^Uj) . The corresponding 
Hamiltonian H gets the form 

H =*% (_ ! J) - i(o'(<) +fJ'(f)) (I J) + 

Once again we emphasize that although the Hilbert 
space of the composite system is infinite-dimensional, the 
prescribed time evolution belongs to the state space V 
over the two-dimensional invariant subspace "Hi of the ex- 
citation number operator E = a^a + (l/2)a z + 1/2, where 
only one excitation is present with certainty. In this phys- 
ical context our method generates a Hamiltonian model 
reduced to V. It is easy to check that a reduction of the 
time-independent Hamiltonian model 

ffjc = Hg(a^a- + aa+) (22) 

has a matrix representation coincident with the model 
given by (|2T|) . The Hamiltonian iJjc represents the 
well known Jaynes-Cummings model of a two-level atom 
resonantly coupled to a single electromagnetic field 
mode, taken in the interaction picture with respect to 
Ho = huj (a) a + \<Jt) ■ Even in this example, the same 
mathematical construction would be valid if we substi- 
tute |0,+)(0,+| and |1,-)(1,-| in (l9j with the two ba- 
sis projectors |+)(+| and |— )(— | of a single qubit. In this 
case ([21]) would describe the interaction of the qubit with 
a particular magnetic field. 

V. CONCLUSION AND DISCUSSION 

In the paper we have proposed a method of recon- 
structing a Hamiltonian from a prescribed time evolution 
of a quantum state p. We suppose that the quantum 
system is governed by the von Neumann equation, and 
so every evolution is represented by a unitary operator. 
This implies that the action of the group of evolution op- 
erators over the state space is equivalent to the adjoint 
action of the unitary group over Hcrmitian matrices. As 
a result the conditions Tr p k (t) = const for all fc^ 1, de- 
termining an orbit of the adjoint action of the unitary 
group, are fulfilled. For every such state p it is possible 



to construct the class of evolution operators realizing the 
prescribed evolution, and this evolution is located within 
a single orbit. 

In the theory of open quantum systems one deals with 
reduced density matrices. They are governed by more 
complicate equations than the von Neumann one. These 
equations, known as quantum master equations, are dif- 
ferent for different systems. Evidently, the dynamics of 
an open quantum system, in general, can not be repre- 
sented by a unitary operator. Nevertheless, we hope that 
the inverse dynamical problem for reduced density matri- 
ces may be solved within the group-theoretical approach 
reported in this paper. 
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